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Abstract 

Under assumption of singular behavior of invariant charge as{q^) 
at ~ and of large behavior, corresponding to the perturbation 
theory up to four loops, a procedure is considered of smooth match- 
ing the /3-function at a boundary of perturbative and nonperturbative 
regions. The procedure results in a model for as for all > with di- 
mensionless parameters being fixed and dimensional parameters being 
expressed in terms of only one quantity Aqcd- The gluon condensate 
which is defined by the nonperturbative part of the invariant charge is 
calculated for two variants of "true perturbative" invariant charge, cor- 
responding to freezing option and to analytic one in nonperturbative 
region. Dimensional parameters are fixed by varying normalization 
condition asim^) — 0.29, 0.30, 0.36. It is obtained that on the 
boundary of perturbative region uq = as((?o) — 0.44, the procedure 
results in nonperturbative Coulomb component acouiomb — 0.25, the 
nonperturbative region scale go — 1 GeV, the model parameter a ~ 
(0.42 GeV)^ which suits as string tension parameter, the gluon con- 
densate appears to be close for two variants considered, K ~ (0.33 - 
0.36 GeV)4 (for ^^(m^) = 0.33). 
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1 Introduction 



A consistent description of interaction of fundamental fields both at large 
and short distances is one of the most important problems of QCD. The 
strength of the interaction is defined by invariant charge as^q"^) (the running 
coupling constant), which satisfies the renormalization group equation. The 
purpose of the present work consists in a formulation of a model for as(g^) 
for all > 0, which is appropriate for description both of perturbative and 
of nonperturbative phenomena and needs minimal number of parameters.^ 

We assume that there exists some value go, which characterizes the non- 
perturbative effects scale and the corresponding value ao = ctsiq^) such, 
that for q'^ > q^ and thus as{q^) < Oq the finite loop perturbation theory 
is applicable and sufficient, while nonperturbative effects prevails in region 

< ^Q, which contains also nonphysical singularities of the perturbation 
theory and so here this theory essentially needs an extension of definition. 
In general, the invariant charge and the beta-function as well are depending 
on a renormalization scheme [21 El • For definiteness while performing calcu- 
lations at > we use MS scheme. It is well-known, that the /3-function 
in the perturbative QCD is of the form 

/3pert(/i) = -boh^ - 2hh' - ^h' - hh' + 0(/i6), h = ^4^. (1) 

z 47r 

For = 3 we have values of coefficients 6o = 9, 6i = 32, 62 — 1287.67, 63 ~ 
12090.38 (coefficients bo, bi do not depend on renormalization scheme while 
values 62, ^3 correspond to a choice of MS'-scheme). Expressions obtained 
by solution of Gell-Mann - Low equation 

with the use of are widely used for sufficiently large momenta transfer, 
however they can not be applied in the infrared region. 

As a matter of fact the behavior of at small momenta till now is 
an open question. Lattice methods and SD equations give no an ultimate 
answer. For behavior of the invariant charge as{q^) at — > a number of 
variants are considered (see, e.g. Ref. [3j). Now the most popular variants 
for ^ behavior are: as — > 0, — > 0(1), is strongly enhanced. We 
consider the last possibility, in particular the well-known singular infrared 

^This paper is an extended version of Ref. pp. 
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asymptotic behavior 

-7 T' 9^0' 3 

47r 

for a review see, e.g. Ref. 5 and more recent papers Refs. [H]. Such behav- 
ior of the invariant charge ay in so-called V-scheme corresponds to a linear 
confining quark-antiquark static potential at that g^M"^ = Gvrcr, where a 
is the string tension. Results of some works [7j on the lattice study of the 
three-gluon vertex demonstrate a necessity of taking into account of non- 
perturbative contributions to the running coupling constant being of the 
form Q. In the framework of the continuous QFT additional arguments 
in favor of behavior (jSJ are also presented in paper [S]. Highly similar sin- 
gular infrared behavior has recently developed model for the QCD analytic 
invariant charge 

Asymptotic behavior © occurs provided /3(/i) —h for h ^ oo. We 
consider a possibility of behavior (jS)) and assume the following form of the 
infrared /3-function 

P[h) = -h + z, h> ho, (4) 

where z is a constant and /iq corresponds to the boundary between pertur- 
bative and nonperturbative regions. For h < ho we shall use /3-function (^Q) 
with a finite number of loops taken into account. Our recipe for construction 
of /3-function for all /i > consists in a smooth matching of expressions (0} 
and at point /i = /iq in approximations of the perturbation theory up to 
four loops. The demand of the /3- function and its derivative to be continuous 
uniquely fix free parameters z and ho of the global /^-function (the matched 
one). Note, that the presence of parameter z in Eq. Q which corresponds 
to Coulomb contribution in invariant charge gives a possibility of smooth 
matching. 

We are willing to build the model for the invariant charge, which precisely 
coincides the perturbation theory in the perturbative region > q^, while 
in the nonperturbative infrared region it provides simple description of main 
nonperturbative parameters. The work is organized in the following way. In 
Section 2 we obtain the matched solutions for the cases of 1 — 4 loops. 
Dimensionless parameters of the model are uniquely defined and further 
the solutions are normalized at the scale of the r-lepton mass, that leads 
to definite values of dimensional parameters A, qo, cr. In Section 3 the 
gluon condensate is calculated. In doing this we consider first a possibility 
of freezing of the perturbative component of Og in the infrared region and 
secondly a possibility of analytic behaviour of this component in the infrared 
region. Section 4 contains concluding remarks. 
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2 Construction of One — Four-Loop Matched So- 
lutions for All > 

For an illustration let us consider the most simple one-loop case. Conditions 
of matching give two equations 

-bohl = -ho + z, 



-2boho = -1. (5) 

The solution of set @ reads 

We shall normalize perturbative solution 

bolnx A^^^ 

by value 47r/io, that gives 

xo = Qo/^QCD = e^ (8) 

where e = 2.71828... . Imposing on as{q'^) the natural condition to be contin- 
uous at q'^ = Qq, we may normalize nonperturbative solution of equation @ 

as{q^) = Att ( ^ + z] , < ql (9) 



by AirhQ as well. As a result we obtain 

C ='^ = ql{hQ- z), CO = C/Aqcd = Mho - z). (10) 

Equations (|l()j) are correct for one — four loops, for one- loop case cq = 
e^/46o- For final fixation of the solution for all > we need to define 
^QCD by normalizing the solution, say, at point q'^ = m^, where nir = is 
the mass of the r-lepton. 

For one-loop case we have following simple formulae for quantities under 
consideration 

"0 = 27r/6o, acouiomb = 7r/6o, 
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qo = cAqcd, o- 



3e' 



-A 



fill 



The gluon condensate K for frozen perturbative constituent which is given 
here for completeness wih be calculated below. For nj = 3 one has oq = 
0.698, acouiomb = 0.349, go = 2.72Aqcd, ct = (IMAqcd?, K = 
(0.980Aqc£))'*. For normalization condition, e.g. as(m.^) = 0.32 one obtains 
^QCD = 0.201 GeV and we have only qualitative agreement for quantities 
considered. As we shall see further on, rather poor results in the one-loop 
case are marginal, while results for many-loop cases turn to be much more 
promising. 

Let us consider multi-loop cases. Solution h{q'^) of equation © for L = 
ln(g^/A^) — > oo reads as follows 

46? 



boL 



1 r ^ 



In^L-lnL - 1 + 



bob2 
8bl 



In'^L 



+ - 



1 



b^nb 



+ 



2 - 
1 

l4 



36o&2 
"86f 



InL 



(12) 



2 166f 

Keeping in the expression terms with powers of logarithms in denominators 
up to the first, the second, the third and the fourth, we fix the 1 — 4-loop 
approximations of the perturbation theory for running coupling constant. It 
may be written in the form 



Q,(g2) = 4^/i(g2) 



47r 

— a(x), 

00 



1 



a{x) = 

mx 



In(lnx) 



+ b' 



In^ X 

3^1"-^ 5 1^2 (In x) 



In^(lnx) In(lnx) 



In'^x 



In X 



+ 



In"* X 



In^(lnx) 



, Infln x) 
+ (3k + 1)-V^ + 



Here x = j 



In^ X 2 In^ X In^ x 

and coefficient are defined as follows 



In X 



(13) 



26i 

-1 + 



K 



bob2 
8bj ' 

'0^3 



bU: 



166f 



(14) 
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Coefficients b, k, R depend on n/.With = 3 we have b ~ 0.7901, n ~ 
0.4147, K ~ —1.3679. In the case of the two-loop approximation for pertur- 
bative ag we have the following set of equations 



bohl + 2bihl = ho 



1 



Inxo 



The set for three loops reads 



2boho + Gbihl 
ln(ln Xq) 

In^ Xq 



1, 

: boho, 



Co , 
\- z 

Xq 



(15) 



bohl + 2bihl + 



2bQho + 6bihl + 2b2hl 



1, 



1 Jn(lnxo) 2 
m xo In^ Xq 



The set for four loops reads 



In^(lnxo) ln(lna;o) 



In^ Xq 
Co 



+ 



In Xq In Xq 



bQho, 



Xq 



+ Z = hn 



(16) 



1 



In Xq 

In^(lnxo) 
In^ Xq 



bQhl + 2bihl + -^hl + 63/1^ = h 

2bQho + 661 /ig + 262 /iQ + 563/10 
^ ln(lnxo) ^ ^2 



1, 



In^(lnxo) In(lnxo) 



In^ Xq 
5 In^(lnxo) 



In*^ xo 



In Xq 
+ (3k + 1) 



In^ Xq 

In (In xo) 
In^ Xq 



+ 



K 



+ 



In Xo 
R 



In^ Xo 



'0"0; 



Co , 
\- z 

Xq 



hn 



(17) 



From sets of equations (|T^ - H17() we find values of Hq, z, xq and cq. They 
are presented in Table ^ Taking into account the existing data ^1 ITTl 
I12| we finally fix the momentum dependence of solutions for a number of 
values of the running coupling constant at r-lepton mass scale nir with the 
effective number of flavors Uf = 3. The values of h.QCD corresponding to 
these normalization conditions are presented in Table |2J values of boundary 
momentum qq = ^/xqAqcd are presented in Table El the string tension 
parameter a = (87rco/3)AQ(-,^ is given in Table 01 
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Table 1: The dimensionless parameters h^, z, xq = Qq/Aq^j-), cq = C/Aq 
on the number of loops, Uf = 3. 





1-loop 


2-loop 


3-loop 


4-loop 


ho 


0.0556 


0.0392 


0.0356 


0.0337 


z 


0.0278 


0.0215 


0.0203 


0.0197 


Xo 


7.3891 


7.7763 


10.2622 


12.4305 


Co 


0.2053 


0.1374 


0.1572 


0.1741 



Table 2: Values of the parameter Aqcd (GeV) on loop numbers and nor- 
malization conditions. Normalization conditions: as(m^) = 0.29, 0.30, 
0.36 with rur = 1.77703 GeV, n/ = 3. 



as{mi) 1-loop 2-loop 3-loop 4-loop 



0.29 


0.1600 


0.3168 


0.2873 


0.2837 


0.30 


0.1734 


0.3370 


0.3069 


0.3026 


0.31 


0.1869 


0.3568 


0.3263 


0.3212 


0.32 


0.2005 


0.3762 


0.3454 


0.3394 


0.33 


0.2143 


0.3951 


0.3642 


0.3573 


0.34 


0.2280 


0.4136 


0.3827 


0.3749 


0.35 


0.2418 


0.4315 


0.4007 


0.3920 


0.36 


0.2556 


0.4490 


0.4184 


0.4087 



3 Gluon Condensate 

Let us turn to calculation of the gluon condensate. Its value is defined by 
the nonperturbative part of Os- We have (see, e.g. the third of Refs. 0) 

oo 

K ^< as/TT : G^,G^, :>= J dq' q^a^^t{q^). (18) 



We define the nonperturbative part of the invariant charge as a difference 
of total invariant charge and its perturbative part (to all orders). In our 
approach the nonperturbative contributions are presented for < Qq only. 
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Table 3: Values of the parameter qq (GeV) on loop numbers and normal- 
ization conditions. Normalization conditions are the same as in Table [21 



as(m^) 1-loop 2-loop 3-loop 4-loop 



0.29 


0.4350 


0.8833 


0.9203 


1.0002 


0.30 


0.4713 


0.9397 


0.9831 


1.0667 


0.31 


0.5081 


0.9949 


1.0452 


1.1323 


0.32 


0.5451 


1.0490 


1.1065 


1.1967 


0.33 


0.5824 


1.1018 


1.1667 


1.2598 


0.34 


0.6198 


1.1533 


1.2258 


1.3216 


0.35 


0.6572 


1.2034 


1.2838 


1.3820 


0.36 


0.6947 


1.2522 


1.3405 


1.4409 



So we have 



^=-3/ dq^1^Mq^)-'^ptiq^)) = ^ j dq^ q^{as{q^) - a^t{q^)). (19) 


Let us consider two variants of the invariant charge perturbative component 
behaviour. 

3.1 Freezing of perturbative component 

For the beginning we define the perturbative part in nonperturbative region 
basing on the assumption of freezing of apt at small q^ (see, e.g. Ref. 1131). 
That is we assume 

apt(g^) = as{ql) = 47r/io, q^ < ql- (20) 
Using expressions ©, (f^ . we have 



= -{ho - z)xlAl^CD- (21) 

Expression H21|) is valid for each of 1 — 4-loop approximations of in 
perturbative region and ratio Kf^/ Kqqj^ does not depend on normalization 
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Table 4: String tension parameter \fa (GeV) on loop numbers and normal- 
ization conditions. Normalization conditions are the same as in Table [21 



as(mi^) 1-loop 2-loop 3-loop 4-loop 



0.29 


0.2098 


0.3398 


0.3297 


0.3426 


0.30 


0.2274 


0.3615 


0.3522 


0.3654 


0.31 


0.2451 


0.3827 


0.3745 


0.3878 


0.32 


0.2630 


0.4035 


0.3964 


0.4099 


0.33 


0.2809 


0.4239 


0.4180 


0.4315 


0.34 


0.2990 


0.4437 


0.4392 


0.4527 


0.35 


0.3170 


0.4629 


0.4599 


0.4733 


0.36 


0.3351 


0.4817 


0.4802 


0.4935 



conditions for Og. The values of K^l /Kqcd calculated with the aid of 
expression (j21jl are presented in Table 13 For normalized solutions values of 

Table 5: Dimensionless quantities K^J^ / Aqcd, Uo, / Kqcd on loop 

numbers, nj = 3. 





l-loop 


2-loop 


3-loop 


4-loop 


yo 

e 

kU^/^qcd 


0.9799 
1.0 

7.3893 
0.9373 


0.8977 
0.8299 
9.3702 

0.8494 


0.9952 
2.2691 
4.5225 
0.9353 


1.0709 
2.8710 
4.3297 
1.0025 



the gluon condensate K^J are presented in Table IHl It is seen that the two 
— four-loop results are quite stable with respect to loop numbers and for 
normalization condition as(m^) = 0.33 one has Kfj- = (0.355 - 0.383 GeV)^, 
which is close to the conventional value jT^ of the gluon condensate (0.33 
GeV)l 

3.2 Analytic behavior of perturbative component 

Let us consider another variant of definition of the perturbative part of for 
< Qq- Namely instead of freezing ()20|) we shall assume "forced" analytic 
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Table 6: Values of the gluon condensate K^/^ (GeV) on loop numbers 
and normalization conditions. Normalization conditions are the same as in 
Table El 



as{m'^) 1-loop 2-loop 3-loop 4-loop 



0.29 


0.1568 


0.2844 


0.2859 


0.3038 


0.30 


0.1699 


0.3025 


0.3054 


0.3240 


0.31 


0.1832 


0.3203 


0.3247 


0.3439 


0.32 


0.1965 


0.3377 


0.3437 


0.3635 


0.33 


0.2099 


0.3547 


0.3624 


0.3827 


0.34 


0.2234 


0.3713 


0.3808 


0.4014 


0.35 


0.2369 


0.3874 


0.3988 


0.4198 


0.36 


0.2504 


0.4031 


0.4164 


0.4377 



behavior of apt in this region, 

apt(g') = aan('?'), < ql (22) 

The main ideas of the analytic approach in quantum field theory, which 
allows one to overcome difficulties, connected with nonphysical singularities 
in perturbative expressions, were proposed in Refs. |151 116| . The analytic 
approach is successfully applied to QCD [TJj. The forced analytic running 
coupling constant is defined by the spectral representation 

oo 

aan(2/) = - / — : — pier), (23) 


where spectral density p{a) = 9oan(— cr — ^0) = 9a(— a — iO). For the 
perturbative solutions a{x) of the form (|13j) the two-loop analytic running 
coupling constant and its nonperturbative part were studied in Ref. |18j . 
the three-loop case and the four-loop case where studied in Refs. [TP and 
Refs. |20j, respectively. Let us write the spectral density up to the four-loop 
case. 

P^^^(-) = ^' (24) 
pi^){a) = p('Ha) - yr^^ [2irtFr{t) - - vr^) F,{t)] , (25) 
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p(3)(cT) = p(2)(cT) + ^^^^ [vT {3t' - n') {Flit) - Fiit)) 
-2t {t^ - Stt^) Fi{t)F2{t) - TT {3t^ - vr^) Fi{t) +t [t^ - 3tt^) F2{t) 

+TTK{3t^ -tt"^)] , (26) 



63 



+47rt {t^ - vr^) (Ff (t) - 3Fi(t)F|(t)) - lOvrt (t^ _ (^2^^) _ ^2^^)^ 
+5 (t^ - 6^2i2 ^ ^4-) F^(^t)F2{t) + 47r (1 + 3k) t (t^ - Fi(t) 

- (1 + 3k) (t^ - eTT^t^ ^ ^4^ ^ 4^-^ ^^2 _ ^2-)j ^ ^27) 

Here t = ln{a), 

Fi(t) = iln(t2 + vr2), F2(t) = arccos ^=L=. (28) 
Solving the equation^ 

aa.n{y) = boho, (29) 

we find values yo for Oan, being defined by formulas H23() - (PS)) with the 
use of values /iq obtained above. Further we find dimensionless quantity 
^ = (Aan/Agc"!))^ = xq/uq. Values of yo and in dependence on the num- 
ber of loops are presented in Table El In Table values of Aan = Qo/y/Vo 
are presented in dependence on the number of loops and on normalization 
conditions at = m^. 

Let us turn to the gluon condensate. In the considered method of defi- 
nition of the perturbative part of ag in the nonperturbative region we have 

go io 
-^an = ^ j dq^ <f{a^s{q^) - aan(g^)) = ^ J dq'^ <f + z 





^While solving this equation it is convenient to use the method of Refs. fl9l . in which 
Ci^{y) is represented as a series in inverse powers of y. 
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Taking into account the spectral representation (|23|) and performing inte- 
gration in Eq. (|3U() . we obtain 



12 



an ^2 V""^ ' 2' 



12A|n 



da p{a) 



yo - cr In ( 1 + 



yo 

a 



(31) 



Substitution a = exp(t) leads to the following representation of the gluon 
condensate 



Kan — J^QCD 



12 

vr 



2 V ^0 - 77 j 



I2e 



dtp{t) {yoe* 



-e2Mn(l + yoe-*)}]. (32) 

The expression in the square brackets of (|32|) does not depend on values of 
as(m^), values of the ratio kH'^ /Aqcd, which are obtained by numerical 
integration in formula (|32() , are given in Table \S\ In Table [7| and Table |S1 
values of the parameter Aan and of the gluon condensate kH'^ are presented, 
respectively. 



Table 7: Dependence of the parameter Aan (GeV) on loop numbers and 
normalization conditions. Normalization conditions are the same as in Ta- 
ble El 





1-loop 


2-loop 


3-loop 


4-loop 


0.29 


0.4350 


0.9696 


0.6109 


0.5903 


0.30 


0.4714 


1.0315 


0.6526 


0.6296 


0.31 


0.5081 


1.0921 


0.6939 


0.6682 


0.32 


0.5451 


1.1515 


0.7345 


0.7063 


0.33 


0.5824 


1.2094 


0.7745 


0.7435 


0.34 


0.6198 


1.2659 


0.8138 


0.7800 


0.35 


0.6572 


1.3210 


0.8522 


0.8156 


0.36 


0.6947 


1.3745 


0.8899 


0.8504 



For cases from one loop up to four loops the behavior of the running 
coupling constant 0^(0'^) for all > is shown in Fig. ^ Here the behavior 
of the analytic coupling constant aan((/^) for < q^, which defines the 
perturbative part of as{q^) in this region, is also shown. In Fig.[21in addition 
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Table 8: Values of the gluon condensate ii'an (GeV) in dependence on loop 
numbers and normalization conditions. Normalization conditions are the 
same as in Table |2I 



as{m'i) 1-loop 2-loop 3-loop 4-loop 



0.29 


0.1500 


0.2691 


0.2687 


0.2844 


0.30 


0.1625 


0.2862 


0.2870 


0.3033 


0.31 


0.1752 


0.3031 


0.3052 


0.3219 


0.32 


0.1880 


0.3195 


0.3230 


0.3403 


0.33 


0.2008 


0.3356 


0.3406 


0.3582 


0.34 


0.2137 


0.3513 


0.3579 


0.3758 


0.35 


0.2266 


0.3666 


0.3748 


0.3929 


0.36 


0.2395 


0.3814 


0.3914 


0.4097 



to as{q'^) its nonperturbative part OnptC?^) is shown, which turns to be zero 
for > ^Q. Normalization condition for all curves in Fig. ^ Fig. [21 is 
as{ml) = 0.32. 

4 Conclusion 

Starting from the well-known perturbative expression for the /3-function 
for small values of the coupling constant h and the behavior @ for large 
values of the coupling constant, which corresponds to the linear confine- 
ment, we have constructed the model /3-function for all h > 0. To control 
dependence of the results on the number of loops we simultaneously consider 
cases corresponding to perturbative /3-function for 1 — 4 loops. 

While constructing the matched /3-function we assume it to be smooth at 
the matching point, that leads to the invariant charge being smooth together 
with its derivative for all > 0. The normalization of the invariant charge, 
e.g. at rrir fix it thoroughly. Then value as(m^) ~ 0.33 corresponds to value 
of the parameter o"^/^ ~ 0.42 GeV of our model which fits the string tension 
parameter of the string model 21 . 

The obtained invariant charge is applied to study of the important phys- 
ical quantity, the gluon condensate. In doing this we consider two variants of 
extracting of the nonperturbative contributions from the overall expression 
for the invariant charge. The first variant assumes "freezing" of perturba- 
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tive part of the charge in the nonperturbative region (f' < q^,^ while for the 
second one we choose the analytic behavior of the perturbative part in this 
region. As we see from Tables IHl El the first variant leads to values of the 
gluon condensate being somewhat larger, than that for the second variant. 

Emphasize, that for ^^(m^) = 0.33 values of the gluon condensate for 
both variants are quite satisfactory. Namely the value for the second variant 
practically coincides the conventional value jl4j . while for the first variant 
it is only slightly higher, K^^'^ = K^/^ ~ 0.36 GeV. Note, that other im- 
portant parameter, the nonperturbative scale qo for the same normalization 
condition also turns to be of a reasonable magnitude, go — 1-17 GeV (see Ta- 
ble OJ. We may conclude, that the present model consistently describes the 
most important nonperturbative parameters. From this point of view the re- 
sults support a possibility of singular infrared behaviour of the invariant 
charge which contains both perturbative and nonperturbative contributions. 
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Figure 1: Running coupling constant as{q'^) for all > and analytic 
coupling constant aa,n{Q^) for < qo (the corresponding curves are the 
lower ones). Normalization conditions: agirn^) = 0.32, Q;an(9o) ~ (^siQo)- 
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